function [PosFinale, ValFinale] = pendulum (APAcc,VerAcc,Fs,xIn,xFin)

%% Da Rafael et al 2007

if (xIn*128-2*Fs)<0
    Inizio = 0;
else
    Inizio = xIn*128-2*Fs;
end
if (xFin*128+2*Fs)>length(APAcc)
    Fine = length(APAcc);
else
    Fine = xFin*128+2*Fs;
end
       
APAcc = APAcc(xIn*128-2*Fs:xFin*128+2*Fs);
%Filtro a 2.5Hz per estrarre la componente principale dall'accelerazione
%anteroposteriore. Nel nostro caso CH3
[b,a] = butter(4,(10/Fs),'low');
Armo = filtfilt(b,a,APAcc);

%Trovo i massimi di questa accelerazione IC
[maxpos,maxval] = findMax(Fs,Armo,10,Fs*2,(xFin-xIn)*Fs+Fs*2,0,1);

figure;
plot(Armo);
hold on
plot(maxpos,maxval,'rx');

%Trovo i minimi locali nell'intorno dei massimi, sull'accelerazione
%anteroposteriore non filtrata.

PosFinale = 0;
ValFinale = 0;

